Novel hyperbolic metamaterials based on multilayer graphene structures 
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We suggest a new class of hyperbolic metamaterials for THz frequencies based on multilayer 
graphene structures. We calculate the dielectric permittivity tensor of the effective nonlocal medium 
with a periodic stack of graphene layers and demonstrate that tuning from elliptic to hyperbolic 
dispersion can be achieved with an external gate voltage. We reveal that such graphene structures 
can demonstrate a giant Purcell effect that can be used for boosting the THz emission in semicon- 
ductor devices. Tunability of these structures can be enhanced further with an external magnetic 
field which leads to the unconventional hybridization of the TE and TM polarized waves. 

PACS numbers: 42.70.-a, 42.79.-e, 78.67.Bf, 73.20.Mf 
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A hyperbolic medium is a special class of indefinite 
media [l| described by the diagonal permittivity tensor 
with the principal components being of the opposite signs 
which results in a hyperbolic shape of the isofrequency 
contours 0, Such media have a number of unique 
properties including negative refraction 0, 0] and sub- 
wavelength imaging One of the possible realizations 
of hyperbolic media is a periodic metal-dielectric nanos- 
tructured metamaterial where the hyperbolic nature of 
the isofrequency curves appears due to the excitation 
of the near-field plasmon Bloch waves [H, 0]- Hyper- 
bolic metamaterials have been realized for optical, in- 
frared, and microwave frequency ranges. Realization of 
the THz hyperbolic media could allow to boost otherwise 
slow THz radiative transitions in semiconductor devices 
which would lead to the development of a new class of 
THz sources. 

Graphene, a two-dimensional lattice of carbon atoms, 
exhibits a wide range of unique properties [a-llCf. Sur- 
face plasmons excited in individual graphene sheets have 
been extensively studied, both theoretically [lll - [l7j and 
experimentally 3 III- 

In this Letter, we suggest a novel class of hyperbolic 
metamaterials where individual graphene sheets are sep- 
arated by host dielectric slabs, as shown schematically in 
Fig. [TJ It is easy to notice an analogy between a graphene 
sheet placed inside a dielectric medium and a thin metal 
waveguide imbedded into a dielectric matrix, which also 
supports localized surface plasmon polaritons. Assuming 
this analogy, we may expect that a periodic lattice of the 
graphene sheets may behave like an effective hyperbolic 
medium due to the coupling between the surface plas- 
mons localized at the individual graphene sheets (20j ] . Im- 
portantly, surface plasmons in graphene have low losses 
and strong localization in the THz region. Indeed, as 
we demonstrate below, a periodic structure of graphene 
layers creates a novel type of metamaterial with strong 
nonlocal response and hyperbolic properties of its disper- 
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FIG. 1. (Color online) Geometry of the effective hyperbolic 
metamaterial composed of a periodic structure of graphene 
layers with conductivity a separated by dielectric slabs with 
the thickness d and dielectric constant e. 



sion curves for TM-polarized waves in the THz frequency 
range and superior characteristics such as a giant Purcell 
effect and tunability by a gate voltage or magnetic field. 

It is important to mention that the periodic layered 
structure shown in Fig. Q] resembles a natural graphite 
which is known to exhibit many interesting properties in 
the ultraviolet frequency range 21]. However, in the case 
of graphite, the 7r-orbitals of the carbon atoms in the in- 
dividual graphene planes strongly overlap which results 
in the modification of the electron band structure and 
large nonradiative losses. In the metamaterials we sug- 
gest here the electronic structure of the graphene sheets 
remains unperturbed. 

We consider the structure shown in Fig. [T] composed of 
graphene layers separated by dielectric slabs. Rece ntly 
a graphene based hyperlens has been proposed in [22] . 
However, the structure considered in this paper resem- 
bles the wire medium metamaterial and our structure is 
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similar rather to the mutlilayered meta-dielectric meta- 
materials. First, we employ the transfer matrix approach 
and analyze the dispersion properties of electromagnetic 
waves for both polarizations. We describe a graphene 
sheet by macroscopic parameters as a conducting sur- 
face defined by the frequency dependent conductivity 
<j(lu) |2j|. Electromagnetic boundary conditions at the 
graphene plane can be written as a matrix M connecting 
the tangential components of the electromagnetic field 
from both sides, 
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so that the transfers matrix for the structure period can 
be defined as T — SMS, where the matrix S is the trans- 
fer matrix of a dielectric layer. Then, we obtain the dis- 
persion relation for the Bloch waves in the metamaterial 
structure in the form cos(Kd) = (Tn + T22V2, i.e. 



TE : cos(Kd) 
TM : cos(Kd) = cos(k z d) 



cos(k z d) - sin(fc z d), 
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sin(fc z d), 



(2) 
(3) 



where k z = (ek^ — /3 2 ) 1 / 2 , j3 is the wavevector in the 
plane, K is the Bloch wavenumber, and fco = lj/c. 

As is clear from Eqs.([2]) and ©, the dispersion prop- 
erties of the structure are defined by the graphene con- 
ductivity <j{uj) which depends on the chemical potential 
fi. Thus, changing /1 with the external gate voltage we 
can tune the isofrequcncy contours of the metamaterial. 
For example, for the chemical potential the isofrequency 
curves of TM polarized waves are elliptic whereas they 
become hyperbolic for the gate voltage to 10 mV (see 
Fig. [2]). The transition from the elliptical to hyperbolic 
regime has been studied recently in metamaterials [3], 
and the graphene-based metamaterials offer a simple way 
for the realization of such transitions. 



As has been established recently [25[ , dispersion prop- 
erties of nanostructured metal-dielectric metamaterials 
may demonstrate strong optical nonlocality due to ex- 
citation of surface plasmon polaritons, and these effects 
are more pronounced when the thicknesses of metal and 
dielectric layers are dissimilar. Accordingly, we expect 
graphene metamaterials to exhibit strong nonlocal prop- 
erties having one of the layers vanishingly thin. 

To derive the effective parameters characterizing the 
graphene metamaterial, we employ the nonlocal homog- 
enization procedure [26j j . Within this approach, we as- 
sume that the structure is excited with a harmonic exter- 
nal current J ~ exp[i(/3a; + Kz — cot)], and we calculate 
the electrical displacement averaged over the period of 
the structure. After lengthy calculation employing the 
techniques from Ref. [26j], we obtain the nonlocal dielec- 
tric permittivity tensor e n i oc , 
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■ TM, V G =0 
TMV G =10mV 
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FIG. 2. (Color online) Isofrequency contours of the graphene 
metamaterial for the different values of the gate voltage and 
a fixed frequency of 4 meV. 



with all nondiagonal elements vanishing and the remain- 
ing three components 




e xx = E-2a\l-2a^f{p,K) 
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where a = 27rlm(cr) / (kodc) , and the function /(/3, K) has 
the form, 



f(P,K) 



k,Qdsm(k z d) 
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2k z [cos(Kd) - cos(k z d)} (k 2 z - K 2 ) ' 



(8) 



For small values of j3 and K, the function f(/3, K) 
vanishes, and the components of the effective dielec- 
tric permittivity tensor become wave independent as in 
the case of quasistatic approximation and local media, 
^xx ^yy ^ 2ct. 

Figure Eta) shows the dependence of the dielectric ten- 
sor component e xx = e — 2a. on the in-plane wavevector 
P for a fixed frequency of 40 meV, whereas the predic- 
tion of the local approach is shown with a straight dashed 
line. The component e yy deviates slightly from its local 
value e — 2a, and it is not shown in the figure. 

We observe that the nonlocal dielectric permittivity 
becomes negative for small values of pd. Because the 
component e zz is always positive, this results in a hy- 
perbolic shape of the isofrequency contours. However, 
for larger values of P the local approximation becomes 
invalid, and the component e xx changes its sign to pos- 
itive, so an effective hyperbolic medium is transformed 
into elliptic one. This is a direct manifestation of strong 
nonlocal properties of the graphene metamaterials. 
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FIG. 3. (Color online) (a) Dependence of the component e xx 
on the in-plane wavevector /3. Frequency is 40 meV; local 
approximation result is shown with a straight dashed line. 
Insert shows the dependence on K. (b) Frequency dependence 
of the component e xx for different values of j3. 



Figure EJb) shows the dependence of the component 
e xx vs. frequency. A solid line corresponds to the case 
of (3 — 0, and it coincides with the result of the local 
approximation. In this case, e xx may take very large 
negative values. However, there exist again strong dis- 
crepancies between local and nonlocal theories for the 
case of nonzero /?. We observe that in this case there is 
a certain frequency corresponding to the pole of expres- 
sion ([7]) for e xx below which the dielectric permittivity 
becomes positive. 

To illustrate some consequences of these results, we 
calculate the Purcell factor of a point dipole placed on 
top of the graphene metamaterial which characterizes the 
enhancement of the radiation efficiency as compared to 
the free space. In our structure we expect large Purcell 
factors for low frequencies. The expressions for the Pur- 



cell factors for both dipole orientations have the form Q , 
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R \\ = 4 Re / £3 r TM, (10) 
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where signs || and _1_ correspond to the parallel and per- 
pendicular orientations to the interface, respectively, and 
tt e , ttm are the reflection coefficients from the semi- 
infinite layered structure. The Purcell factor is shown 
in Fig. [4] We observe that the low frequency range is 
characterized by large Purcell factors for TM-polarized 
waves. This agrees with our observation of hyperbolic 
isofrequency contours for the low frequency range in TM 
polarization. We also notice that there are no large Pur- 
cell factors values for the TE polarization in low-loss fre- 
quency range. Moreover, for the large frequencies the 
large values of Purcell factors can be observed for both 
polarizations. This increase in Purcell factor is connected 
with a step- like increase of the non-radiative losses in the 
graphene (proportional to the real part of conductivity) . 
Here we should mention the effect of the temperature on 
the spontaneous emission enhancement. The tempera- 
ture defines the steepness of the step-like increase of the 
real part of the conductivity. In order to minimize the 
non-radiative losses in the hyperbolic media region we 
should imply the condition /j/^bT) « 1. In the THz 
frequency range, this condition corresponds to the tem- 
peratures of about 10 Kelvin. 
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FIG. 4. (Color online) Purcell factor vs frequency for both 
TE and TM polarizations. Temperature is 4 K. 

We have also obtained analytical expressions for the 
Ff^ in two limiting cases. When the conductivity of 
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graphene is small, the Purcell factor takes the form, 
3ire 2 ( c \ 3 ( eckod 



R 



TM 
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27r|Im(cr)| 



exp 



27r|Im(<T)| 



> (12) 



We notice that the Purcell factor is approaching zero ex- 
ponentially when conductivity vanishes. When the con- 
ductivity is large, the expression for R?[ reads 
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TM 
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so that in the limit of infinite conductivity the Purcell 
factor vanishes. 

Thus, the Purcell factor vanishes very large and very 
small conductivities. Thus, we should observe some max- 
imum at the intermediate case. We have revealed that 
the maximum value of the Purcell factors is achieved for 
a = e/2. We notice that a — e/2 corresponds to epsilon- 
near-zero metamaterial. In this regime, the media is 
characterized by almost plane isofrequency contours and 
all evanescent modes can be excited. 

In the low frequency region the losses are small for suf- 
ficiently low temperatures. Thus, no contribution to the 
Purcell factor in the low frequency region is due to the 
losses. From the other hand, in our calculations we did 
not account for nonlocal effects in graphene conductiv- 
ity which should become significant for the in-plane vec- 
tors of the order of the Fermi wavevector: kp — [i/vf, 
where vf is the Fermi velocity of graphene electrons, and 
fi is the Fermi energy. We have checked how the non- 
locality would modify the obtained results by perform- 
ing integration in Eqs. (|9]).(fl0|). and (fTT|) to k F - We fix 
the frequency to 44 meV and change the period of the 
structure. We find that for small periods the effective 
model sufficiently overestimates the values of the Pur- 
cell factor. However, when the period is increased the 
values of the Purcell factors for both models are almost 
equivalent. We have also derived an analytical condition 
of good agreement between local model and the model 
which performs integration only over a finite range of 
wavevectors: coth(fid/(2vF )) < 2Tr\lm(cr)\fi/ (ujvf)- A 
good agreement between local and nonlocal models al- 
lows to expect strong THz radiation enhancement in mul- 
tilayered graphene metamaterials. 

We also study the properties of this novel metamate- 
rial in the presence of a strong constant magnetic field. 
In this case, the surface conductivity of graphene is de- 
scribed by a 2 x 2 tensor with nondiagonal component 
equal to the Hall graphene conductivity <jh- Further- 
more, in the presence of strong magnetic field the lon- 
gitudinal and Hall parts of the conductivities are gov- 
erned by the electron transitions between the discrete 
Landau levels 
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In addition, for the case of a single 
graphene sheet the magnetic field couples the TE and TM 
eigenmodes resulting in the emergence of new magneto- 
plasmonic eigenmodes [28[. 



For the layered structure, we observe the similar cou- 
pling between the TE and TM modes. Indeed, the Bloch 
vectors of the system obey the following equation, 



or * A + B HA -Bf ^alsm 2 {k z d) 
cos{K lt2 d) = — — ± \ . 



. 2ma kg . , n 

A = cos(fc z a) — sm(fc z d), 

c k z 

B = cos(k z d) —^sm(k z d). 

c k e 



(14) 



When the coupling term ~ er# vanishes, Eq. (fT4|) simpli- 
fies to Eqs. @ and ([3]). In the presence of the magnetic 
field, the TE and TM polarized modes become coupled 
which results in complicated shapes of the isofrequency 
contours. To illustrate these results, we plot the eigen- 
modes isofrequency contours in the presence of constant 
magnetic field of 1 T applied along the z axis. 
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FIG. 5. (Color online) (a) Frequency dependence of the Hall 
and longitudinal conductivities in the presence of static mag- 
netic field of 0.1 T. Temperature is 4 K. (b) Isofrequency con- 
tours for the frequency 100 meV for the first (red solid line) 
and second (blue dashed line) hybrid electromagnetic modes. 

Figure [5ja) shows the dispersion of the longitudinal 
and Hall conductivities in the vicinity of the resonance 
corresponding to the transition between the 0-th and 1- 
st Landau levels. Figure [5jb) shows the isofrequency 
contours of the electromagnetic eigenmodes at the fre- 
quency 100 meV. We observe that, although the imag- 
inary part of the longitudinal conductivity is negative 
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and we would expect the elliptic contours for both the 
eigenmodes, the coupling via the Hall conductivity term 
changes the properties dramatically, and the isofrequency 
contours have a complicated shape, which is neither el- 
liptical nor hyperbolic. The strong dependence of the 
isofrequency contours on the magnetic field suggests a 
new degree of freedom for engineering optical properties 
of these novel metamaterials. 



Finally, we comment briefly on feasibility of the pre- 
dicted effects. It is technologically achievable to make 
layered structures consisting of several graphene sheets 
separated by dielectric layers. For fabrication of single, 
double and multiple graphene layers, or two-dimensional 
layers of other materials, the most popular approaches 
are cleavage technique [29( as well as CVD [30| and 
MBE [3l[ methods. There exist several methods of layer 
transfers, such as wet and dry transfer pro cedures with 
PMMA, PMGI or water-soluble layer [33-[3i]. These 
methods remain a lot of residue at the later after transfer, 
but the surface can be cleaned with Ar/H2 annealing. In 
a number of studies [32U34j , heterostructures composed 
of several graphene layers, boron nitride, or molybde- 
num flakes have been already demonstrated. Hexagonal 
Boron Nitride (h-BN), conventionally used as a host me- 
dia, can be employed as dielectric with the permittivity 
4 and higher in the THz frequency range. 



In conclusion, we have demonstrated that a multilayer 
periodic structure composed of graphene layers can oper- 
ate as a novel type of hyperbolic metamaterial in the THz 
frequency range, where tuning from elliptic to hyperbolic 
media can be achieved by applying an external gate volt- 
age. We have predicted a giant Pureed effect for a point 
dipole source placed on top of such a metamaterial be- 
ing achieved through coupling of effective plasmon modes 
with infinite wavenumbers. We have also demonstrated 
a novel possibility of tuning the dispersion properties of 
the metamaterials with external magnetic field and the 
emergence of a family of novel hybrid eigenmodes in the 
graphene metamaterials. 
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